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Abstract 

Lusignoli and Vignati have recently pointed out that it is in principle possible 
to directly detect the cosmic antineutrino background (CIaB) by using the rather 
stable isotope 163 Ho as a target, which can decay into 163 Dy via electron capture 
(EC) with a very small energy release. In this paper we calculate the rate of the 
relic antineutrino capture on 163 Ho nuclei against the corresponding EC decay rate 
by taking account of different neutrino mass hierarchies and reasonable values of 
# 13 . We show that such flavor effects are appreciable and even important in some 
cases, and stress that a calorimetric measurement of the CuB might be feasible in 
the far future. 
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The standard cosmology predicts that relic neutrinos and antineutrinos of the Big 
Bang must survive today and form a cosmic background, but their average temperature 
is so low (T v ~ 1.945 K) that their existence has not directly been verified [Tj. Among 
several proposals for a direct detection of the cosmic neutrino background (Cz/B) [2], the 
most promising one should be the neutrino capture experiment by means of radioactive 
/3-decaying nuclei [3] [H]. A direct measurement of the cosmic antineutrino background 
(CuB) could similarly be done by using radioactive nuclei which can decay via electron 
capture (EC) [T2], but it seems much less hopeful. 

Lusignoli and Vignati have recently pointed out an interesting and noteworthy way 
to directly detect the CuB [13J. It takes advantage of the rather stable isotope 163 Ho as 
a target, which undergoes an EC decay into 163 Dy with a small energy release (Q ~ 2.5 
keV). The signature of this relic antineutrino capture process is located on the right- 
hand side of the spectral endpoint of the EC decay, and their interval is detectable if the 
experiment has a sufficiently good energy resolution. In practice, it is promising to do a 
high-statistics calorimetric experiment to measure the 163 Ho spectrum, probe the absolute 
neutrino mass scale [13] and even detect the Cz/B. 

In this paper we shall calculate the rate of the relic antineutrino capture on 163 Ho 
nuclei against the corresponding EC decay rate by taking account of different neutrino 
mass hierarchies and reasonable values of the smallest neutrino mixing angle 8 13 . Such 
flavor effects, which were not considered in Ref. [13], are found to be appreciable and 
even important in some cases. We shall first look at the energy spectrum of the EC 
decay of 163 Ho by examining how sensitive the ratios of the peak heights and the fine 
structure near the spectral endpoint are to the neutrino mass hierarchy, and then present 
a detailed analysis of flavor effects on the detection of the Cz7B by using a target made of 
this isotope. Although our numerical results are mainly for the purpose of illustration, we 
stress that such a direct laboratory measurement of the CuB deserves further attention 
and investigations and might be feasible in the long term. 

The key point of the relic antineutrino capture on 163 Ho nuclei is to capture the 
relic electron antineutrinos in a way without any energy threshold. In the three-flavor 
scheme the |z/ e ) and |z/ e ) states are superpositions of three neutrino and antineutrino mass 
eigenstates, respectively: 

We) = Y, V ekWk) , 
k 

We) =E V ekWk) , (1) 
k 

where V ek (for k = 1,2,3) denotes an element in the first row of the 3x3 neutrino mixing 
matrix V. We have \V el \ = cos9 12 cos9 13 , \V e2 \ = sin 9 l2 cos 9 l3 and \V e3 \ = sin 6^3 in the 
standard parametrization of V [15] , where 9 l2 — 34° and 9 l3 < 12° as indicated by a 
global analysis of current neutrino oscillation data [16]. Of course, 9 13 = 0° and 9 13 7^ 0° 
may have very different consequences in neutrino phenomenology, as we shall see later on. 
Now let us consider the EC decay of the isotope 163 Ho [T7] : 

163 Ho + e^ shcll) -> 163 Dy* + v e -> 163 Dy + E i + u ei (2) 

where e^ shell ) denotes an orbital electron from the i-th shell of holmium-163, and E i is 
the corresponding binding energy of the electron hole in dysprosium- 163. If the Q-value 
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of this EC decay is defined as the mass difference between the two atoms in their ground 
states, then the energy spectrum of the outgoing neutrinos will be given by a series of 
lines at Q — E v Note that Q has not well been determined and its value is possible to 
range from 2.3 keV to 2.8 keV [13]. In what follows we shall take Q = 2.5 keV as a 
default value unless otherwise explicitly mentioned. Taking account of the Breit-Wigner 
resonance form of the atomic levels [17], one may express the energy spectrum of the 
detected EC events as follows: 
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EC 'I* (Q-T c )J2\K k \y(Q-T c f-mle(Q-T c -m k ) 



dT c 4tt 2 

T 1 

x E n, <7< A 2 B t ^ • ^ _ Ei)2 - pp , (3) 

where A EC denotes the rate of this EC decay, T c is the so-called "calorimetric" energy 
which measures the difference between Q and the neutrino energy [T5] , Gp = G F cos 9 C 
with 6 C ~ 13° being the Cabibbo angle of quark flavor mixing, m k is the mass of V k (for 
k — 1, 2, 3), and the theta function Q(Q — T c — m k ) has been imposed on Eq. (3) to ensure 
the kinematic requirement. In addition, n i denotes the fraction of occupancy of the z-th 
electron shell, C i is the nuclear shape factor, fi i represents the Coulomb amplitude of the 
electron radial wave function, and B i is an atomic correction for the electron exchange 
and overlap. The Q- value of this EC decay is so small that only those electrons from M 1 , 
M 2 , N 1 , N 2 , 1 , 2 and P 1 levels can be captured [IB] . In an excellent narrow- width 
approximation [18] the integral of Eq. (3) yields 
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V = -A E E »i ^ (Q - E t ) \V ek \^(Q - E t y - m\ , (4) 

i k 

where % runs over all the possible electron shells of 163 Ho that involve its EC decay, and 
k runs over three neutrino mass eigenstates. 

Given a positive Q-value, a thresholdless capture of the incoming relic antineutrinos 
on the EC-decaying 163 Ho nuclei may happen through 

V e + 163 Ho + e^ shcll) 163 Dy* 163 Dy + E % . (5) 

The capture rate for each mass eigenstate V k (for k = 1, 2, 3) hidden in the flavor eigenstate 
V e can be given by 

\ = ^ % fc I V ek 1 2 E^Q A 2 B t ■ {Evk+ Q_ 1 Ei) 2 + T 2 /4 . (6) 

where % and £L are the number density and energy of u k , respectively. The standard 

k k 

Big Bang cosmology predicts n v = n-^ ~ 56 cm -3 today for each species of relic neutrinos 
and antineutrinos P, but possible gravitational clustering of massive v k and v k around 
the Earth could make their number densities much bigger [19]. For simplicity, we shall 
not consider the gravitational clustering effect in the subsequent analysis. Note that the 
de-excitation energy of unstable 163 Dy* in Eq. (5) is in principle monoenergetic for each 
mass eigenstate V k (i.e., T k = E v + Q). Convoluted with a finite energy resolution in 
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practice, the ideally discrete energy lines of the final states in Eq. (5) must spread and 
then form a continuous spectrum. As usual, we consider a Gaussian energy resolution 
function defined by 

R ( T , T k) = J- 



2na 



(T - T k f 



2a 2 



(7) 



where T is the overall energy of an event detected in the experiment. Using A to denote 
the experimental energy resolution (i.e., the full width at half maximum of a Gaussian 
energy resolution for the detected events), we have A = 2y/2 In 2 o w 2.35482 a. Then the 
differential antineutrino capture rate reads 
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In the mean time, the energy spectrum of the EC decay should also be convoluted with 
the Gaussian energy resolution. So we rewrite Eq. (3) as 
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The above two equations allow us to calculate the rate of the relic antineutrino capture 
on the EC-decaying 163 Ho nuclei against the corresponding background (i.e., the EC de- 
cay itself) in the presence of flavor effects, which are characterized by both the neutrino 
masses m k and the neutrino mixing matrix elements \V ek \ (for k — 1, 2, 3). 



We proceed to evaluate the flavor effects on both the relic antineutrino capture rate 
and the corresponding EC decay rate with the help of Eqs. (8) and (9). In our numerical 
calculations we take Am^ m 7.6 x 10~ 5 eV 2 and (Am^j m 2.4 x 10" 3 eV 2 together with 
9 12 ~ 34° as typical inputs [16]. The impact of 9 13 on A F and A EC can be examined by 
allowing its value to vary within 0° < 9 l3 < 12° [IB] . Depending on the sign of Am| 1; two 
mass patterns of three active neutrinos are possible: one is the normal mass ordering with 
m 1 < m 2 = \J m\ + Am^ < = \ fm\ + |Am§]J ; and the other is the inverted mass 

ordering with m 3 < mi = \jm\ + |Am§j] < m 2 = yfm% + |Am§]] + Am^ . In either 
case the absolute neutrino mass scale (m l or m 3 ) is unknown, but its upper bound is 
expected to be of 0(0.1) eV up to 0(1) eV as constrained by current experimental and 
cosmological data [151 120]. As for the values of those parameters relevant to the atomic 
levels of 163 Dy, we refer the reader to Ref. [13] and references therein. The energy levels 
of the captured electrons are assumed to be fully occupied (i.e., n i = 1) and their binding 
energies and widths can be found in Table 1 of Ref. [33]. The atomic corrections for the 
electron exchange and overlap are neglected (i.e., Bi ~ 1), and the ratios of the squared 
wave functions at the origin (i.e., Pf / can be found in Table 2 of Ref. [13]. Because 
the nuclear shape factors C { are approximately identical in an allowed transition [17], they 
can be factored out from the sum in Eqs. (3), (4), (6), (8) and (9). Finally, one should 
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note that the numerical results of A F and A EC can be properly normalized by means of 
the half-life of 163 Ho through the relation A EC T-,y 2 = In 2, where T x , 2 ~ 4570 yr. Then the 
distributions of the number of capture events and the number of background events are 
expressed, respectively, as 

diV s 1 dA F In 2 _ 



dT A EC dT T 1/2 
diV B 1 dA EC In 2 



dT A EC dT T 1/2 



N T t (10) 



for a given target factor N T (i.e., the number of 163 Ho atoms of the target) and for a given 
exposure time t in the experiment. 

Let us first look at the energy spectrum of the EC decay of 163 Ho. Our main concern is 
the flavor effects on the ratios of the peak heights and the fine structure near the endpoint 
of the energy spectrum. The numerical results are shown in Fig. 1, where 9 13 = 10° has 
typically been input. Fig. 1(a) and Fig. 1(b) show that two ratios of the peak heights 
are completely insensitive to the lightest neutrino mass m , which is identified with m 1 
for the normal mass hierarchy (i.e., m 1 < m 2 < m 3 ) or with m 3 for the inverted mass 
hierarchy (i.e., m 3 < m 1 < m 2 ). Given m < 0.5 eV, the actual variations of two ratios 
are only at the C(10~ 6 ) level. Hence it is hopeless to obtain any information about the 
absolute mass scale of three neutrinos or their mass hierarchy by measuring these ratios. 
However, it seems very promising to acquire knowledge of the Q-value from these ratios, 
as one can clearly see in Fig. 1(c) and Fig. 1(d). In Fig. 1(e) and Fig. 1(f) we illustrate 
how the ratios of the peak heights change with the energy resolution A for given values 
of m and Q. We see that the changes are at the percent level if A < 1 eV is taken. 

The fine structure near the spectral endpoint of the EC decay of 163 Ho is illustrated 
in Fig. 2, where 6 13 = 10° has also been input. The location of the endpoint is essentially 
determined by the smallest neutrino mass. It might be possible to distinguish between 
the normal and inverted neutrino mass hierarchies, since they have different effects on 
the energy spectrum near its endpoint. Especially for the inverted mass hierarchy, an 
obvious distortion occurs at the position T c ~ Q — m l with m l = \Jm 3 + |Ara§ij due t° 
the fact that the lightest mass eigenstate v 3 hidden in v e is associated with the smallest 
neutrino mixing matrix element \V e3 \ 2 = sin 2 # 13 . In the region of T c > Q — m 1 , only 
the v 3 component contributes to the spectrum. In comparison, the contributions of u l 
and v 2 components to the spectrum are dominant in the region of T c < Q — m 1 . One 
might therefore be able to determine the magnitude of | V^ 3 1 2 if the endpoint distortion in 
the energy spectrum could really be observed. A similar distortion cannot appear in the 
case of the normal neutrino mass hierarchy, simply because the contribution of the lightest 
mass eigenstate u x to the spectrum is of the same order of magnitude as those of v 2 and v 3 . 

Let us now focus on the relic antineutrino capture on 163 Ho nuclei against the EC decay 
background by taking account of the flavor effects and examining the Q- value dependence. 
In our calculations we assume 30 kg 163 Ho as a reference isotope source. The numerical 
results are presented in Fig. 3, Fig. 4 and Fig. 5, where the relic antineutrino capture rate 
is illustrated as a function of the overall energy release T in a possible experiment. Some 
comments and discussions are in order. 



5 



(1) Given 9 13 = 10° as a typical input, Fig. 3 shows the effects of the neutrino mass 
hierarchy on the CvB signature and the corresponding background. The value of the 
finite energy resolution A is chosen in such a way that only a single peak of the signature 
can be seen; namely, at least one peak of the signature is not allowed to fall into the 
shade of the background. Eq. (8) tells us that the contribution of each antineutrino mass 
eigenstate V k (for k = 1,2,3) to the capture rate is located at T = T k = E v ^ + Q and 
weighted by the relic antineutrino number density n v ~ 56 cm" 3 and the flavor mixing 
matrix element |V^ fc | 2 . On the other hand, Eq. (9) tells us that the energy spectrum of 
the EC decay near its endpoint is dominated by the lightest neutrino mass eigenstate 
hidden in V e and sensitive to the energy resolution A. As the smallest neutrino mass (m l 
in the left panel of Fig. 3 or m 3 in the right panel of Fig. 3) increases from to 0.1 eV, 
the CvB signature moves towards the larger T — Q region. In comparison, the shift of 
the corresponding background is less obvious because the smallest neutrino mass and A 
have the opposite effects on the location of the spectral endpoint of the EC decay. Hence 
the distance between the peak of the signature and the background becomes larger for a 
larger value of the smallest neutrino mass, and accordingly the required energy resolution 
A becomes less stringent. Comparing between Fig. 3(c) and Fig. 3(d), for example, one 
can also see that it is more or less easier to detect the CvB in the Am^ < case, where 
the signature is separated more obviously from the background. The reason is simply that 
the dominant antineutrino mass eigenstates V 1 and V 2 have slightly larger eigenvalues in 
this case than in the Am^ > case. Note that more than one peak may emerge in the 
distribution of the relic antineutrino capture rate, as shown in Fig. 3(a) and Fig. 3(b), 
provided the energy resolution is sufficiently good. But the requirement A < m k /2 [6] is 
in general needed so as to resolve the corresponding peak. 

(2) The effect of the unknown neutrino mixing angle 9 l3 on the CiaB detection is 
examined in Fig. 4, where both the Am^ > case with m 1 = and the Am| x < case 
with m 3 = are taken into account. The energy resolution A is chosen in such a way 
that the signature in Fig. 4(b) with m 3 = and 9 13 = 0° should be clearly seen. In the 
Amjj > case, Fig. 4(a), Fig. 4(c) and Fig. 4(e) show that a change of 6*13 in its allowed 
region (i.e., 9 13 < 12°) almost has no effect on both the signature and the background. 
In the Am 31 < case, however, the background turns out to be sensitive to a change of 
9 13 as one can see from Fig. 4(b), Fig. 4(d) and Fig. 4(f). As 9 13 becomes larger, the 
signature essentially keeps intact but the background apparently moves towards the larger 
T — Q region, making the signature partly fall into the shade. Note that the distribution 
of the CvB capture rate may have a slight distortion only when 9 l3 is sufficiently large, 
as shown in Fig. 4(e) and Fig. 4(f). Hence it is in practice almost impossible to probe 
the magnitude of 9 13 by detecting the CuB. 

(3) The Q-value may sensitively affect the absolute capture rate of relic antineutrinos, 
but it does not influence the relative locations of the signature and its corresponding 
background. This property is illustrated in Fig. 5, where the smallest neutrino mass (i.e., 
m 1 in the Am 31 > case or m 3 in the Am 31 < case) has been taken to be zero for 
simplicity. Allowing Q to vary from 2.3 keV to 2.8 keV, we find that the absolute capture 
rate undergoes a decrease of as much as one order of magnitude. The reason for this 
feature can be explained with the help of Eq. (8). The binding energies E i of the seven 
available atomic levels are all smaller than 2.1 keV, and thus the individual capture rates 
at T k = E v +Q are strongly suppressed by the corresponding Breit-Wigner distributions. 



6 



A larger Q-value implies a stronger suppression of the overall capture rate A F , leading to 
a worse signature of the Cz/B in the experiment. That is why the isotope of a smaller 
Q- value is eagerly wanted [32], in order to make this interesting Cz/B detection method 
more feasible in the foreseeable future. 

In summary, we have examined the flavor effects on a possible detection of the Cz/B 
by using a target made of 163 Ho nuclei. On the one hand, we have discussed the en- 
ergy spectrum of the EC decay of 163 Ho by looking at how sensitive the ratios of the 
peak heights and the fine structure near the spectral endpoint are to the neutrino mass 
hierarchy. On the other hand, we have presented a detailed analysis of the effects of 
different neutrino masses and reasonable values of the smallest neutrino mixing angle 8 13 
on the relic antineutrino capture rate and the corresponding EC decay background. Our 
results demonstrate that such flavor effects, which were not considered in Ref. [13J, can 
be appreciable and even important in some cases. 

At least three factors may apparently affect the observability of relic antineutrinos 
in such a calorimetric experiment. The first one is the number of the target particles 
N T , which is the only adjustable parameter for a given measurement to increase the 
total capture events. In our numerical analysis we have taken 30 kg 163 Ho as a reference 
isotope source, just for the purpose of illustration. A much larger target is certainly 
welcome, but it will certainly be a bigger challenge to a realistic experiment. The second 
factor is the finite energy resolution A, which has a crucial impact on the signature-to- 
background ratio. As the precision of A is gradually improved, it should be possible 
to establish a signature of relic antineutrinos just beyond the spectral endpoint of the 
EC decay of 163 Ho in the future. We admit that A = 0.015 eV taken in our numerical 
examples is too optimistic, because it is about two orders of magnitude better than the 
present achievement. If the energy resolution were really perfect, we would even be able 
to resolve the multiple peaks coming from different neutrino mass eigenstates and then 
obtain some information on neutrino masses and flavor mixing angles. But it will be 
a great success even if only a single peak of the signature can in practice be resolved. 
The third factor is the number density of relic antineutrinos n v ^ (for k = 1,2,3) around 
the Earth. Its value might be much larger than n v = ~ 56 cm~ 3 predicted by the 
standard Big Bang cosmology, if the gravitational clustering effect is significant [19] . A 
larger value of n v will always be a good news for the detection of the Cz/B. 

k 

To conclude, the Cz/B and Cz/B can in principle be detected by means of the /3-decaying 
nuclei (e.g., 3 H) and the EC-decaying nuclei (e.g., 163 Ho), respectively. Although the 
present experimental techniques are unable to lead us to a guaranteed measurement of 
relic neutrinos and antineutrinos in the near future, we might be able to make a success 
of this great exploration in the long term. Of course, it would be much harder to identify 
flavor effects in the signatures of relic neutrinos and antineutrinos. Let us reiterate that 
a direct laboratory measurement of the Cz/B and Cz/B is extremely important and thus 
deserves more attention and our best efforts. At least from a historical point of view, 
nothing is completely hopeless in neutrino physics. 

This work was supported in part by the China Postdoctoral Science Foundation under 
grant No. 20100480025 (Y.F.L.) and in part by the National Natural Science Foundation 
of China under grant No. 10875131 (Z.Z.X.). 
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Figure 1: The energy spectrum of the EC decay of 163 Ho: the ratios of the peak heights 
M 1 /M 2 and M 1 /N 1 changing with the lightest neutrino mass m (equal to m 1 in the 
m| x > case with a solid curve or to m 3 in the m% l < case with a dashed curve), the 
Q-value and the energy resolution A. Here 9 13 = 10° has typically been input. 
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Figure 2: The energy spectrum of the EC decay of 163 Ho: the fine structure near its 
endpoint in the mfj > (top panel) or < (bottom panel) case. Here Q = 2.5 keV 
and 8 13 = 10° have typically been input. 
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Figure 3: Effects of neutrino masses on the relic antineutrino capture rate as a function of 
the overall energy release T in the Am 31 > (left panel) or Am| x < (right panel) case. 
The solid and dashed curves represent the CuB signature and its background, respectively. 
The value of the finite energy resolution A is chosen in such a way that only a single peak 
of the signature can be seen. Here 9 l3 = 10° has typically been input. 
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Figure 4: Effects of 6 l3 on the relic antineutrino capture rate as a function of the overall 
energy release T in the Am^ > case with m 1 = (left panel) or in the Am^ < case 
with m 3 = (right panel). The solid and dashed curves represent the CuB signature and 
its background, respectively. The value of the finite energy resolution A is chosen in such 
a way that the signature in Fig. 4(b) with #13 = 0° should be clearly seen. 
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Figure 5: Effects of the Q-value on the relic antineutrino capture rate as a function of the 
overall energy release T in the Am^ > case with m 1 = (left panel) or in the Am| x < 
case with m 3 = (right panel). The solid and dashed curves represent the Cl^B signature 
and its background, respectively. The value of the finite energy resolution A is chosen in 
such a way that only a single peak of the signature can be seen. Here Q = 2.3 keV, 2.5 
keV and 2.8 keV together with 8 13 = 10° have typically been input. 
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